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1. INTRODUCTION AND SUMMARY 

Work Assignment #2 to Contract Number DAADO5-72-C-0284 is concerned with 
development of a Glint Model for the Cobra AH-1G to analyze glint off the 
windows, Sun glints off the rotor hub, fuselage and windows in level flight 
and while banking allow the visual observer to detect low flying aircraft at 
much longer ranges than is possible in nonglint conditions. Reflections and 
glints from the rotor blades, rotor hub amd fuselage have been adequately 
taken care of by use of a non-reflecting light-dispersing paint. This task 
is then concerned with window areas. 

The objectives of the task are to develop a simulation model of the Co- 
bra AH-1G for use in predicting the visual detectivity of the aircraft by a 
ground observer due to sun glints from the windows; to exercise this model to 
determine probabilities of detection vs. angles of incidence of sunlight; and 
to use the model to determine the best sunshade configuration. These objec- 
tives were met by the model described in Section 3 which used the window con- 
figuration described in Section 2. The model included the parameters of sun 
angle, sun intensity, atmospheric transmission, sun-aircraft relationship, 
detection sensitivity of the human eye, background intensity and contrast 
between background and aircraft (to determine detection ranges for observer 
position relative to the aircraft). A 50% probability of detection criterion 
was used to determine the ranges. The required reflectivities for the win- 
dows were computed as a function of angle of incidence for reducing detection 
range to the nonglint detection range of the helicopter and for reducing de- 
tection range to 1500 meters. 

The results showed that to reduce detection range to the nonglint detec- 


tion range of the helicopter, the required reflectivity varied from less than 
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.1% to more than 30% over all angles of incidence. To reduce detection range 
to 1500 meters, the required reflectivity varied from less than .001% to more 
than .5% over all angles. The most stringent requirements were due to large 
glints off the relatively flat overhead window. If this window was removed 
from the model, the low end requirements were relaxed by an order of magnitude 
to less than 2% and less than .01% respectively. These data indicate that 
large glints occur at all angles of incidence and that a combination of 
improved window design, improved window coatings and a sunshade system 
may be required to reduce glints to an acceptable level. 

It is recommended that additional work in this area be performed 
on shutter designs and then use the model developed under this task as a 
basic tool in evaluating the shutter redesigns. The model should be modified 
for future use to consider cluttered and uncluttered backgrounds with high 
and low background reflectance, This modification is a minor one. In addi- 
tion, the criteria for use in evaluating shutter designs must be determined 
and agreed upon. 

The model as developed to this point may also be used in evaluation of 
other helicopter windows by inputting the appropriate window information into 
the model, 
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2. THE HELICOPTER WINDOWS 

The Cobra AH-1G has five windows, two on each side and an overhead. The 
location of the windows is described by five outline drawings (listed below) 
and their contour is described by Bell Helicopter drawing 209-B3 Revision A. 


Window Outline Drawing Contour Drawing 
Pilot's Door 209-030-516 

(Pilot's left side) 

Pilot's Window 209-030-507 209-83 
(Right side) 

Gunner's Door 209-030-515 

(Left side) 

Gunner's Window 209-030-508 209-B3 
(Right side) 

Overhead Window 209-030-509 209-53 


The contour drawing describes the mold used in making the windows, Al- 
though it describes only the right hand side and overhead windows, the left 
side windows are the mirror image of the right hand side windows. The format 
of the drawing is a series of twenty-seven full scale curves taken at various 
flight stations from 53.5 to 169.5 along the longitudinal axis of the Cobra 
at intervals of at least .5" and no more than 6", 

In order to describe the windows for use in the model, xyz information 
was generated for 726 points on the right side windows and for 12 points on 
the overhead window using the curves on the Bell drawing. Points were gen- 
erated only for those portions of each curve which fell within the outline of 
the windows as determined by the five outline drawings noted above. 

The center of the coordinate system for the x, y, z values was taken 


as the 0, 0, O point for Flight Station 0, Bulkhead Line 0 and Water Line O. 
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This point is shown in figure 2-1 which was taken from figure 12-2 of the Cobra 


Manual? » This system assigns negative x values to Flight Station numbers . 
positive y values to left Bulkhead Line mmbers and positive z values to 
Water Line numbers, The x y z points for the side windows were determined 
by plotting points at approximately 1" intervals along each of the twenty- 
seven curves of the contour drawing and noting their x, y, z values. Although 
this process could be done by a digitizing plotter with a printout in thou- 
sandths of an inch; it was done by hand for this task due to lead time require-— 
ments for the digitizer. x,y,z values were measured to the nearest hundredth of 
an inch, with an estimated accuracy of +,01 inches. The errors in blueprint 
distortion were estimated to be 1/16" in 36" or less than .5%. The right 
hand window was assigned negative x values, negative y values and positive 
z values. Corresponding points on the left hand window had identical x and z 
values and identical y values with a positive rather than negative sign. 

The overhead window was flat along y and curved in x and z with flats 
at the front and top. It was described by twelve points, with the first and 
second describing the front flat, the eleventh and twelth describing the 
top flat and the second through eleventh describing the curved portion. 


loperator's Manual, Army Model AH-1G Helicopter,T™ 55-1520-221-10, Headquarters, 


Department of the Army, 19 June 1971. 
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Figure 2-1. Helicopter Reference Points 


2-3 
UNCLASSIFIED 


UNCLASSIFIED 


© 


3. THE COMPUTER MODEL 


3.1 General 


The computer model was programmed in Fortran IV and run on the Univac 


1108 digital computer at the Westinghouse and Defense Center. The model output 


includes the following: 


a) 


b) 


c) 


d) 
e) 


f) 


g) 


h) 


i) 
J) 


A listing of 726 side window surfaces and 12 center window 
surfaces with their x,y,z coordinates, x,y,z components of 
their normals and effective area. 


A listing of the 1464 surfaces of both side windows and the 
center window by surface number with x,y,z coordinates for 
each surface, 


An Area Contour Plot showing the projection of all normals in 
zenith and azimuth coordinates. 


A listing of the azimuth and zenith coordinates of each normal. 


A Brightness Plot showing the relative brightness of reflected 
rays in zenith and azimuth coordinates. 


A Helicopter Detection Range contour plot showing the detection 
range of the helicopter along each reflected ray in zenith 

and azimuth coordinates. 

A Window Detection Range contour plot showing the detection range 
of the glint along each reflected ray in zenith and azimuth co- 
ordinates. 


A listing of the azimuth and zenith coordinates, the angle of 
incidence and the glint detection range for each reflected ray. 


A listing of the required reflectivity for each angle of incidence. 


A plot of required reflectivity vs. angle of incidence. 


A flow chart of the model is shown in figure 3-1. Paragraph numbers are 


inserted in the flow chart to show where discussion on the various parts of 


the model can be found, In order to keep machine time and thus computer 


costs down, the model was divided into three segments. The first segment 


computed normals and effective areas. Its output is a) above. This output 
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INPUT XYZ OF INPUT XYZ OF 

‘726 POINTS ON 12 POINTS ON 

SIDE WINDOW OVERHEAD WINDO 
CALCULATE NORMALS CALCULATE NORMALS 
AND EFFECTIVE AREAS OF AND EFFECTIVE AREAS 0 
726 SIDE WINDOW POINTS 12 OVERHEAD WINDOW POINTS 


LIST OF XYZ, NORMAL, AND AREA 
DATA FOR EACH POINT 


DATA CARD FOR DATA CARD FOR 
EACH POINT (726) EACH POINT (12) 


CALCULATE BRIGHTNESS FACTOR 
FOR EACH REFLECTED RAY IN 
« , PB COORDINATES 


CALCULATE AREA IN “, B 
COORDINATES FOR EACH OF 
1464 NORMALS (726+726+12) 


CALCULATE HELICOPTER BRIGHTNESS AREA CONTOUR LISTING OF «,B 
DETECTION RANGE PLOT PLOT FOR EACH NORMAL 
(WITHOUT WINDOW GLINT) 


HELICOPTER CALCULATE WINDOW WINDOW DETECTION RANGE 
DETECTION DETECTION RANGE CONTOUR PLQ 
RANGE 
CALCULATE REQUIRED LISTING OF « ,B, RANGE, AND 
REFLECTIVITY ANGLE OF INCIDENCE FOR EACH 
NORMAL 


REQUIRED REFLECTIVITY 
LISTING AND PLOT 


Figure 3-1. Flow Chart for Cobra Model 
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is fed into the secom segment which computes an Area Contour Plot and an 


associated %,B (zenith, azimuth) plot. Its output is b), c) and d) above. 
This output is required only once for a given window configuration. The out- 
put of the first segment is also fed into the third segment which computes 
brightness factors, detection ranges and required reflectivities for various 
sun positions, visibility conditions and reflectivity requirement criteria. 
Its output is e) through g) above. A listing of the three segments is pro- 
vided in Appendix A. 
3.2 Normal Computation 

The normal computation segment of the model was performed in two dif- 
ferent programs, one for the right side window and one for the overhead win- 
dow. Both programs computed tangents on the window and used these tangents 
to compute x,y,z components of the normals and effective areas. Since this 
data doesn't change for the same window, it was computed only once, and out- 
put on IBM cards which were then used as an integral part of the next program 
segments, The data generated for the right side window applied to the left side 
window when the signs of the y position and the y component of the normal 
were changed from - to +, 

The 726 points on the side window were grouped in 27 sets of fram 1 to 
40 points which described contour curves of equal x values. Figure 3-2 shows 
a representation of several of these lines. Data for each of these points 
was read into the program. The data was then ordered in first x and then z. 
Then tangents along the x curve were computed by computing the vectors between 
adjacent points along the x curve. The tangents were described in x, y and z 
values, Each point, P 


X,Y, 2’ 


the vector from Py )2 to Px y 2-1? 


had tangents Tl and T3 computed where Tl was 
the point with the next lowest z value on 
3-3 
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Figure 3-2, Tangents on Side Windows 
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that x curve, and T3 was the vector from PY, ~ to oa y,2ztl? 


the next highest z value on that x curve. At the end points, artificial vec- 


the point with 


tors of length .001 were introduced to provide continuity. Tangents T2 and 
Th across X curves were computed by finding the closest points on either side 


of the x curves and computing the vectors to them from P On the ex- 


x7 a 
treme points, artificial vectors of length .001 were introduced to provide 
continuity. An example of the four tangent vectors for one point is given in 
figure 3-2, 

To obtain the normals, cross products between Tl and T2, T2 and T3, 

T3 and T4 and T4 and Tl were computed. The result of each cross product 
computation was an x value, y value, z value and an area. The x value of 
the normal was obtained by dividing the sum of the four x values by the sum 
of the four area values. The y and z values were obtained in a similar man- 
ner, The result was a unit normal which was unaffected by the artificial 
tangents since they had a small length and thus a near zero area. A repre- 
sentation of the four subnormals, N12, N23, N34 and N4l and the resultant 
normal N is shown in figure 3-3. 

The angle between each normal and its four subnormals was obtained by 
taking the dot product of the subnormal and the normal. The result was the 
cosine of the angle. This angle was then divided into 1/8° and the result 
multiplied by the area obtained from the cross product computation to obtain 
the effective area for the subnormal. The four subnormal effective areas for 
each area were then added to provide a total effective area for each normal. 
This area was a function of the area of the sector defined by the normal, 
the radius of curvature of the sector and the +1/4° angular subtense of the 


sun. 
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Figure 3-3. Normals on Side Window 
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A similar but simplified method was used for the twelve points on the 


overhead window. Two tangents, Tl amd T3, were computed for each point using 
the adjacent points, A third tangent, T2 of x0, y=16, z=0 was used to rep- 
resent the window width. The normals and effective areas were then generated 
using two cross products, Tl x T2 and T2 x T3. 

The output of the two programs was a listing of 738 points by their 
x,y,Z% position with x,y,z components of their unit vectors and their effec- 
tive areas. This output is presented in Appendix B. An IBM data card was 
generated for each point for use in the next segments of the model. 

3.3 Area Contour Plot 

This segment of the model converts the data of 738 window sectors into 
1464 window sectors covering both side windows (726 sectors each) and the 
overhead window (12 sectors). 

The Area Contour Plot shows the projections of the normals in azimuth 
and zenith coordinates for a trimmed helicopter. The zenith angle of the nor- 
mal is called « and the azimuth angle is called $8 on the plot. Figure 3-4 
shows the ~ , B angular relationships. The zenith and azimuth angles are 
quantized in 3° increments and the sum of all areas within +14° of each 
quantized zenith and azimuth point are plotted on the Area Contour Plot. Num- 
bers on the plot are two place integers representing 10 log, (Area x 1000). 
This representation was chosen to guarantee that all areas could be repre- 
sented by a positive two place number. 

The output of this segment of the model is a listing by surface number 
of the x,y,z, coordinates of each of the 1464 normals, a listing by surface 
number of the & , PB coordinates of each normal, and the Area Contour Plot. 


By assigning surface numbers and listing the ™ , i of each normal, it is 
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Figure 3-4. Alpha, Beta Coordinate System 
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possible to find the point or points on the windows which a given ™% ,P or 


area corresponds to. 

The surface x,y,z listing, the > , F listing and the Area Contour Plot 
are presented in Appendix B., 

The data from this segment of the model is generated only once because 
it does not change for a given window configuration. This information is 
then used in the third segment of the model. 

3.4 Brightness Plot 

This segment of the program computes the reflected sun rays off the 
helicopter window and computes a brightness factor which is used later to 
compute detection range. A Brightness Plot showing relative brightness as 
a function of the azimuth and zenith position of the reflected rays is gen- 
erated, 

Each sun position relative to the helicopter's trim position is consic- 
ered as one case. For each case a new brightness plot is generated showing 
all of the reflected rays. The x,y,z components of the reflected ray are 
computed by equations (1), (2) and (3). 

(1) RX = SX - 2x SDN x XN 

(2) RY = SY - 2x SDN x YN 

(3) RZ=SZ-2x SDNx IN 
where : 

SX = sun rays x component = -cos (sun azimuth) x sin (sun zenith) 


SY = sun rays y component = -sin (sun azimuth) x sin (sun zenith) 


~” 
N 
lI 


sun rays 2 component = -cos (sun zenith) 
XN = x component of normal 
YN = y component of normal 
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ZN = Z component of normal 
SDN = dot product of normal and sun ray = cosine of angle of incidence 
between sun ray and surface normal 
When SDN is negative, the angle of incidence is greater than 90° and the 
normal under consideration is in the shadow of the window. For these normals 
there is no reflected ray and no brightness factor. 
The sun ray vector is defined as a unit vector and the normal vector is 
a unit vector, thus the reflected ray vector is a unit vector. This vector 
is used to calculate the ~ , B position of the reflected ray in a manner 
similar to that of the & , B position of the normal vector in section 3.3, 
with the , B position quantised in 3° sectors. 
The brightness factor is computed by equation (4). 
(4) BF = area x SDN x RO 
where: 


area = the area associated with the normal as computed in the normal 
computation 


RO = coefficient of reflectivity of the window material (acrylic plastic) 
and the angle of incidence. This coefficient is computed in a 
continuous fraction subroutine which uses data supplied by LWL. 

The highest brightness factor for each quantized < , B is printed on 
the Brightness Plot. Numbers on the plot are two place integers representing 
10 logs, (Brightness Factor x 10’), This representation was chosen to guaran- 
tee that all brightness factors could be represented by a positive two place 
number. 

A typical Brightness Plot (sun zenith = 45°) is presented in Appendix B. 
3.5 Helicopter Detection Range 


Helicopter detection range is the 50% probability of detecting the 
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helicopter with the photopic eye and is calculated along each reflected ray. 
It is used in the determination of required reflectivity. 

Two contrast equations are used in the calculations of detection range. 
The first equation (5) describes the contrast of the helicopter against an 


east horizon sky. The second equation (6) describes the contrast required 
in order to see the helicopter. 


(5) CT, = BH x Be “BH x 100 x e (73°9075 x RI/RVIS) 


where : 
BH = horizon brightness = .0009 w/ om”/ steradian 
ROE = reflectivity of helicopter paint = .15 


RT = detection range (50% probability) of helicopter in nautical 
miles 


RVIS = visibility range in nautical miles 


se +9075) _ two percent minimum contrast requirement 


(6) CT, = 1.57 + 36.5 RI“/A 


where: 


A = area in square feet 


The area of the helicopter is calculated as the visible area of the 
helicopter along a reflected ray as calculated by equation (7). 
(7) AP= RX x35 + RY x450 + RZ x 157.5 


where ; 
35 = frontal area in square feet 
450 = side area in square feet 


157.5 = top area in square feet 
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Equations (5) and (6) are transcendental when solving for R. Thus, equa- 


tion (8) and its derivative, equation (9) were used to solve for R. The equa- 
tions approach a solution as F approaches 0. The pair of equations were iter- 
ated with initial R = 0 and -F/F being used as the correction factor, For 
each following iteration, R was incremented by -F/F' until F/F' was less than 
1 miles, 

(8) F = -CT, +1.57 +36 x RI“/AP 


(9) F' = 3,9075 CT,/RVIS + 2 x 36.5 x RT/AP 


The Helicopter Detection Range Contour is obtained by using an a , & 
plot similar to that for the Brightness Plot except that the maximum Helicop- 
ter Range in nautical miles is plotted for each quantized & ,B position. 

A typical Helicopter Detection Range Contour Plot (sun zenith = 45°, sun 
azimuth = 45°) is presented in Appendix B. 
3.6 Window Detection Range 

The window detection range is the range at which there is a 50% proba- 
bility of detecting the window glint with the photopic eye. 

As with helicopter detection range, two contrast equations, (10) and (11) 
respectively, describe the contrast of the sun glint against an eastern hori- 
zon sky and the contrast required to detect the glint, and difference equation 
(12) and its derivative, equation (13), are used to solve equations (10) and 
(32), 


(10) OG, = BS x EYE x RO -BH x 100xe (-3.9075 RG/RVIS) 
BH 


where: 
BS = sun brightness = 2000 w/cm“/steradian 


EYE = photopic response of eye to the sun=.15 
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RG = detection range (50% probability) of window glint in nautical miles 


(11) CG, = 1.57 +36.5 x RG/ABF 


2 


where: ABF = BF/(RO x 144.) = effective area of window sector in 
square feet 


(12) F = -CG, + 1.57 + 36,5 RG?/ABF 
(13) F' = 3,9075 CG,/RVIS + 2 x 36,5 RG*/ABF 


The Window Detection Range Contour is obtained by using an &% , B plot 
similar to that for the Brightness Plot except that the maximm Window 
Detection Range in nautical miles is plotted for each quantized o , B posi- 
tion. A typical Window Detection Range Contour Plot (sun azimuth = 45°, sun 
zenith = 45°) is presented in Appendix B. 

The relationship between brightness factor and window detection range 
can be shown by letting CG, of equation (10) equal CG, of equation (31). 
After some algebraic manipulation, we obtain equation (14). 


(1h) BF= ____ 36.5% RG’ x BHO 
(BS x EYE - BH) x 100 e (72°%075 RG/RVIS) 


Figure 3-5 shows the relationship for visibility ranges of 5, 10 and 20 
nautical miles. 
3.7 Required Reflectivity 
Required reflectivity is defined as that reflectivity which the windows 
would need in order that the detection range (50% probability) of the glint 
would not exceed a specified range, possibly the detection range of the heli- 
copter itself, or a constant range such as 1500 meters. It is generated as 
a function of the angle of incidence of the sun's rays on the window surfaces. 
The required reflectivity is computed by equating CG, of equation (10) 
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and CG, of equation (11) and solving for RO with RT, the desired range, sub- 
stituted for RG. When helicopter detection range is the desired criteria, 
its value as computed earlier is used. When a constant range is specified, 
that range is used for RT. The resulting equation for required reflectivity 
(ROG) is given as equation (14). 


ABF 
BS x EYE x 100 x e(73+9075 RI/RVIS) 


The above calculation is made for each reflected ray. The angle of incidence 
of each ray is computed as arccosine (SDN) and rounded to the nearest degree. 
Then a lowest value of ROG is found for each angle of incidence from 0° to 
90°, Where there are no reflected rays for an angle of incidence, ROG is 
assumed to be an arbitrary number, 10, which is higher than the actual maxi- 
mum reflectivity of which any passive surface can have. 

A listing and a plot of required reflectivity as a function of angle of 
incidence is generated for each case. The plot includes a count of the num 
ber of rays at each incident angle. At the end of a run of multiple cases, 
a summary listing and plot is printed out. An example of the listing and 
plot of required reflectivities is presented in Appendix B. 
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4. EXERCISING THE MODEL 

The first two segments of the model were used one time to compute the 
window normals and the effective area associated with each normal and to gen- 
erate an Area Contour Plot. These results are presented in Appendix B. 

In the third segment of the program, several sets of runs were made for 
various visibility ranges and required reflectivity criteria. Each set of runs 
included 17 sun position cases as shown in Table 4-1. These sun positions 
relative to the helicopter covered the range of possible sun position- 
helicopter relationships. For example, a helicopter banking at 45° against 
a horizon sun might have a zenith = 135°, azimuth = 90° sun position rela- 


tionship. Six sets of data were run as tabulated in Table 4-2. Different 


Table 4-1. Sun Positions Relative to Helicopter 
Azimuth 


ee 2 ee 2 ee ee ee A. ee 


| 


runs at 20 and 5 nautical miles were made to show the effect of visibility 
range on required reflectivity. Different runs using Helicopter Detection 
Range and a constant 1500 meters showed the effect of choosing the different 
criteria for judging window design. The runs without the overhead windows 
were run because the largest glints were found to occur on the overhead win- 


dows and these large glints overpowered the effect of glints off the curved 
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Table 4-2. Parameter Variation in Runs Made 


a : Reflectiv- 


ity Re- 
quirements 


1 Helicopter Detection 5 Included] Per Ta- | Figure 4-1 
Range ble 4-1 

Helicopter Detection 20 Included| Per Ta- | Figure 4-2 
Range ble 4-1 


RS rer rrr gc 


Helicopter Detection 
Range 


1500 meters 


5 Included| Per Ta- | Figure 4-4 
ble 4-1 

1500 meters Included] Per Ta- | Figure 4-5 
ble 4-1 

1500 meters 5 Excluded| Per Ta- | Figure 4-6 
ble 4-1 


side windows. 


The output for each case in the 6 runs included the brightness plot, 
a listing of normal numbers, X , 8 , angle of incidence and glint detection 
range, and a listing and plot of the required reflectivity vs. angle of in- 
cidence, The summary reflectivity listings and plots for cases 1 through 6 
are shown in figures 4-1 through 4-6. The reflectivity listings and plots 


for the seventeen cases in run #1 are included in Appendix C. 
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5. CONCLUSIONS 

From the data obtained in exercising the model, we can obtain several 
conclusions about the nature of the window glint. These conclusions point 
the way to future work required in designing the window. 

5.1 Effect of Visibility Range on Required Reflectivity 

Figures 4-1 and 4-2 show the results of runs #1 and #2 with visibilities 
of 5 and 20 n. miles respectively and reflectivity requirements based on re- 
ducing detection range to that of the non-glint helicopter, This criteria 
produces more stringent reflectivity requirements for lower visibility ranges 
because the non-glint detection range is affected more by the reduced visi- 
bility than the glint detection range. 

However, if the reflectivity requirement is based on reducing detection 
range to a fixed value, as done in runs #4 and #5 (figures 4-4 and 4-5), the 
more stringent reflectivity requirements occur at higher visibility ranges. 
5.2 Effect of Detection Criteria on Required Reflectivity 

A comparison of the results of runs #1 and #4 (figures 4-1 and 4-4) and 
of the results of runs #2 and #5 (figures 4-2 and 4-5), shows a marked dif- 
ference in requirements based on the selection of either reducing detection 
range to the non-glint detection range of the helicopter or reducing detection 
range to a fixed range (1500 meters in runs #4 and #5) as the detection cri- 
teria, This result points out the importance of selecting the criteria to be 
used in evaluating window modifications or redesign early in any future work 
effort to guarantee the proper evaluation of proposed solutions. 

5.3 Effect of Flat vs. Curved Windows 

Early evaluation of the results of exercising the model showed the lar- 

gest glints came from the large flat sections of the overhead and that as sun 
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position relative to the helicopter changed, these glints caused very strin- 
gent reflectivity requirements over all angles of incidence. When the overhead 
and thus the flats were removed from the models, as in cases #3 and #6 (see 
results in figures 4-3 and 4-6), the more stringent reflectivity requirements 
were reduced by an order of magnitude. 

This implies that any window redesign should attempt to eliminate flats 
and should probably attempt to achieve a minimum value for radius of curva- 
ture. 

5.4 Use of Sunshades 

The goal of using sunshades would be to block off certain sections of 
the windows which caused the largest glints or to limit the angles of inci- 
dence over certain portions of the window to ease the reflectivity require- 
ments. It appears that sunshades would be sucessful in achieving this goal, 
but that they alone would not solve the complete problem of reducing sun glints 
because of the large glints remaining over low angles of incidence, A more 
detailed investigation of the effects of sunshades on both reduction of 
glint and reduction of pilot visibility should be made before a final 


recommendation is made, 
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APPENDIX A 
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DIGITAL COMPUTER PROGRAM 
FOR 
COBRA WINDOW MODEL 
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1 “o5euU008 12.580» (7e21Ue e4497Be «79190» © 35130° 
L eobdeuue LeeVO059 78-0409 e5U812r «72746r «40217 
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mbleeb0e 14650098 9e0U0% 226124" «95062» 61372Ur 
moleco05e 1L4eU609 09-960" «2838501 «936369 616825 
-oled>d0r 13.780» 70-920? eS535411e 2.91094, °19897>» 
moleco00 1356420 71.860r «36888r «88016 225344, 
molecoUe 13-0908 726600" «39191* «857759 o2735e0 
=“olec5s00 12472U0% 7303500 e41871+ «626807 «31374» 
“olecdJUr 1222907 74e000e 6460845e «77104 «36124, 
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-59.u00° 14-640» 62-850U» °©18379» ©98208e-.U0003» 
=5)9.e-Uy0r 14.570» b3-75U» e22164» 97165 ©06327> 
“29eUU09 = 1l4e40095 O4.720" 021944" 1.97013" e08706r 
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—-)9-UUU0r &e9107 74+68009 «701867 «37603 53065» 
-59-UQUOU, BeO0U0U5 75e22000% 2.65809" 131531» 
VDAIACAUN) @ YON) ®Z0N) OXNON) © YNON) @ ZNON) 0D ANCN) (N=7O030 720) / 
o7Tevu04e 13688098 6564500 «236135 «9707390009 
-57euUUs 13.830» 054730» e23493r °©96759» 03861? 
“Q7eUUU9 136590 6647007 635589" «50542 11829» 
“D7eUU098 136290" o7e6600% 635883" «89285" e2302be 
“oJ eure 12-9109 06-600" 237852 ©36647> e2625r 
—~oy/7eJuU0P 12.470» 609 .e-5uUUe e4$2944K, 061913» 0531449 
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— a7 euyUr 11-510» 71-10U% 255237» °©65839" «45174 
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*¥oLANK UUUI00 
BLA L UVicvu56 
BLKe Vloz3uU 


LATcHWAL ReFEKCIWCES (BLUCKe WAME) 
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GUU6 NwUUS 
vuU7 NLOLS 
UuLU NiO2» 
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Quuu UUUV2Z02 LYF YQuu1 00u130 22L 0001 000226 230G 
Uuul QUUZS6 SUL OUU1 000247 4OL 0001 000110 65 
Uuul UUUL171 YeL vv001 00U201 93L 0000 R 000254 A 
UUU¥ R 0199040 & UULO I OOU1L4S4 IA 0004 { QOOOUU ITALPH 
UUuLN I OUUEZ47 IC v0LO0 I 000255 ICON 0000 IJ O0Q00U0 ICONV 
UOuu | OUUEHS [OUT OUuU I OOU25<e Ix 0000 I 000253 IXxx 
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VLiWFEINSION LCONV(1U0) 

JIMENSION IA(0l) 

COMMOW/BLK1L/X (738) 6 Y (738) 0 Z(738) © XN(736) © YN( 738) & ZN( 738) 0 AN(738) 
COMMUN/BULARZZ/ LALPH( 1404) » IBETA( 1464) eRT (1464) 6 RG(1464) 9 B( 1464) 
VATA (CICONV(N) epN=le1l1U0)/2HOLs 2HU02]2 2HOSe 2HOY4? QGHOSe cHOO? QHO7 2 2HOB» 
L2rHu9e CHLOe CHIL» 2HL2e CH13e 2tHH14s 2H1LSs 2H1602H1 79 2H16e 2H19e2H20» 
22HeEl se cH22 0 2H232 2H249 cH259 2H260 2H27 9 2H26 0 QH29 0 CHS 
S2ihobedrid20cH339 2HS4e 2H35¢e 2H360 2H37 2 2H38e 2H39e QH4O 
42H41 9c H42!)9 CH43 9 2H449 CHYS? CH4 Or CHY47 & 2H48 2 QHY49e CHOU 
52:4D19 214529 2H559 2HD4e 2CHD5 2? 2H50e 2HD7 2 2HD8e 2H59e QHOU 
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IF (Io JAliv) eNOeIYYIGU TO 25 
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DVe %J K=nt1 
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Cuuu I Q00S02 IC10 0000 I 000014 IDATA 0000 I 014340 IIWN 
GUUL I 014421 IROG 0000 I 011240 ISU 0000 I 014410 J 
GWUvU IT U1l44.a3 JU OUV0 I 014350 JUD 0000 J 014414 JJu 
GCOLuU Lf vl4422 Uk OU00 1 014346 K 0000 J 014412 L 
UUuY I U14503 WN OV00 I 014364 NN 0000 I 014337 NNORT 
Cuuvu K 014403 k 0O0U0 kK 014336 RALIAN 0004 R 010450 RG 
LUULU KR ULY4uU7 RUG 00Uu0 R 014200 ROGL 0004 R 005560 RT 
LuUUU kK UlLYS72 RX 0000 k 014373 RY 0000 R 014374 RZ 
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LUUU K UL4556 SY 0OvUu00 kK 014357 S2 0000 R 014130 S2t 
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COU R UVOJad? ZN 0000 RK 014367 ZNORM™M 0000 R 006350 2R 
de L COGRA wlLNuOw 
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4 UIMENWSTON XR(01464) 6 YR( 1464) 2 ZR (1464) 
ae VIMENSION 19D0(1464) 
Ue OsMENSION SZE (20) eSAZI (20) sROGL (91) 
ie COrmiMGiNW/BLK1/X (738) 6 Y (758) 6 2Z(738) »XN( 738) 6 YN(738) 0 2N(738) 0 AN( 738) 
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RAUL AN=57-29578 
WNORT=,464 
ITineb 
LOUT=v 
BS=2UUUe 
KUcHelS 
bH=.0uU09 
EYt=eld 
KEALD(CaIN® 105) (KC(K) eK=LeNXY) 
EMD LIN® 105) (YC(K) eK=LeNXY) 
105 FORMAT (6F6.3) 
READCLIN® 49) U0 
4Y FURMAT(I3) 
KEAD(LIN® 48) (SZE(UUD) »SAZI (JUD) » JUD=1+ UD) 
48 FOKMA] (2F 8-2) 
l2 FORMAT(5H BETA) 
S56 FURMAI(5X*°41(3HXXX)) 
39 FORMAT (7X¢ LHU0 8X0 2H159 8X2 2H30 2 8X9 2H45e BX Pe QHOUl BX 2 CH75 2 BX 2 2HIU 
Z7KeSHLOSe 7Xe SHL20¢ 7X0 SHL3S5e 7X 3H150° 7Xe 3H165e 7Xe SH180) 
537 FOKMAT (65X + SHALPHA) 
SU1 FORMAT (47X*e 34HHELICOPTER DETECTION RANGE CONTOUR) 
ve FURMAT(43Xe42HSCALE FACTOR: VALUE=RANGE IN NAUTICAL MILES) 
4U1 FORMAT (49X*+3UHWINDOW DETECTION RANGE CONTOUR) 
UY 6UU 1OA=1°91 
RIVTCLIQGA)=106 
oUU CONTINUE 
Dy 8U wj=1e UD 
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VY ée6U ITOA=1991 
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RVY=SY¥m2e*SUN*YNORM 
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AR (NJ =RX 

YRINJIERY 
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SUNN="SDN 

IF (Suv) 101°101°102 
SD=ATAN2 (SORT (1=—SUN*SDN) » SON) 
SD=SU*RADIAN. 
ISu(NISINT(SD+.5) 

CALL CONFRA(NXY*+XC*e#Y¥CeSDeRQ) 
4 (iv) =AREA*SDN*RO 

B(iW)=5(N) #1-0E+07 
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3 iv) =U 8 
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LaLPrig=u 

[dt TAIN) =U 
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CunTINUE 

ALPHA=ATAN2 (SQGRT(1"ZR(N) *®ZR(N)) ®&ZRIN) ) 
ALPHA=ALPHA*RADI AN 
ITALPH(\)=34IN1T (ALPHA/3.+-5) 
bE TA=ATAN2(YR(N) @XRIN)) 

be TA=BbE TA*RADIAN 

1F (BETA) 44243043 

Int Tal,j2=3* INT(BETA/3S=—.5) 

v0 TO 42 
lbtiTal:,J=3S4INT(BETA/3S.+.5) 
CONTINUE 
dF=pbliww)/7(144.41.04+07) 
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K=Ue 
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IF (F/r ltel)4oUr46lre4ol 
Hou R=ReF/F1 | 
GU TU 462 
40] RTUN) =, 
K=Ue 
ele E=ULXR (3 e9U7TD*EK/RVIS) *CG*¥RKU 
FI-E+1. S7tR*R* 353602 5¥*RO/S BF 
Fl=Se9U7TH*e/KVIS+2 6 #RK¥360 D*RO/ SF 
LF (F/R ate lL) 47004710471 
47U K=ker/F 1 
ov TU ule 
471 KoOIN) =k 
IcA=Iwt (Su+1e5) 
IC(1On)=1C(ITOa) +1 
KUg= (1125743605 *R1 (N) *RT(N) *RO/BF) *BH/ (100. ¥*EXP(-5-9075*RT(N) Z 
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ir (RUG=ROGL(10A) ) 26172620262 
cosa ROUOLIIUA) =RKOG 
LF (ROG=RTOT (IVA) )OU1Ls? 2629262 
oU1 KTUT( 120A) =ROG 
coe COlvIIiINUE 
2UL CONTINUE 
whITeE(LOUT»2Z5u) 
<bv FORMAT (iLbil) 
wRITE (CLOUT? 20U2)S5Zt(M) eSAZI(M) 
wUe FURMAT(4&4XelIHSUN ZENITH=0F60e0°04Xe*L2HSUN AZIMUTH="F 6.0) 
WRITE C IOUT? 2u1) 
ea FORMAT (56A° LSHBRIGHTIVESS PLOT) 
wkiTE (IOUT? 2U3) 
<Uo FORMAT (24Ar36HSCALE FACTOR: VALUE=10X LOG BASE 10 OF» 
24z2ri AREA X RLFLECTIVATY X 10 TO THE 7Tid POWER) 
wRaTE( LOUT+?12) 
J=< 
CALL PLOTT(NNORT? J) 
wkiLTe«(LOUT? 39) 
WRITE (L0UT? 37) 
wkaTe( ~OUTe2bu) 
whiTL(1,OUT» 301) 
wxdTE CLOUT? S02) 
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J=3 
CALL PLOTT(NNOKT? J) 
WRITE (TOUT? 3u) 
WRITE CLOUT? 39) 
wkKiTE( [OUT? 37) 
wKiTEL (JOU e250) 
WRiATEC]TOUT* 401) 
ARLTE(] OUT + 3uUe2) 
WkaTe (TOUT 1l2) 
J=h 
CALL PLOTT (NWORT? JU) 
WhiTh ( 1OUT + 38) 
WKITE( {OUT +39) 
WKiTeE ( [OUT+37) 
WRITE CLOUT? 2b0) 
wKITE CLOUT? 5CU) 
DOU FURMA) (29X*4OHREFLECTED RAY ALPHA+e BETAr ANGLE OF INCI» 
cOuribtivCEs LETECTION RANGE LISTING) 
WUM=1NT(NNORT/4) 
J=u 
JO 34U L=leNUM 
J=JUt1 
LF (Jotuel)dwkKITE(IOUT*® 243) 
e4o FURMAIT(1X04(29H WUM ALPHA WETA ANGLE RANGE) ) 
JI=L*4—5 
Juu=JUt3 
wKLTE CLOUT + 341) (LL»+ LALPH(LL) oe LisETA(LL) pISO(LL) pRGI(LL) pL LSUUe UJJ) 
O44 FOKMAT(1X%04(14¢ 3I oer Fue 00 1X) ) 
LF (Jere 40) WRITE (1OUT» 250) 
Te (Jetrue 40) JEU 
w4u CONT LivuUe 
wRiTe (TOUT? 250) 
wit LTE ( fOUT*? 505) 
905 FORMAT (21X+4OHREQUIRED REFLECTIVITY AS A FUNCTION OF Ap 
2l7riNGLE OF IWCLDENCE) 
WraTk(,OUlTe51uU) 
y1lU FOKMAT (1X9 SHANGLE ? 3X0 1HOe 8Xe LHL BXeLH2+BXe LH3SeBXe LH4eBXe LHS kX 
ClHorSArlH7 e+ BXeltHBr axe iH9) 
YU 3S1lU JU=Aal1*9 
JJJ=JU¥*1U79 
JUII=UUIt9 
JJJJIU=UId*10-10 
Wea Tul ,QUUTs S11) JUUde (ROGL ILL) pLLSUGUUe VU) 
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Wika TE 4 yOUT? S06) 

FURMAT(//) 

INLITE CLOUT» SUD) 

by 85bU J=1291 

Iviad(U)=IC(U)/10 
Ie(JJS=aNT (iC (JU) -Lu*ICl10(U) +21) 
CUNTIiWue 
WR4ATEVLOUT® 449) (IC1LO(U) eU=19°91) 
Mate CLOUT® 449) (IC(J) pJ=1e91) 

FORMAT (7X*eY111) 

Wita TLC GOUT ®#? 430) 

FURMAT (SXe ZHhO) 

VO 44u JJ=19°10 

VU 45U J=ir91 

Jh=i 17d 

RUVEKUGL (J) #1 VEFT09 
IRUG=aiNT(100*ALO0G10(hOUG)) 

ak CIKUGe GE e UA¥*1LIU0) LDATA( J) =IBX 

LF CLRUGeLTe JA*1U0) IDATA(JU) =I68L 

COnwT dive 

LF (Jd eEQeL) WRITE (CLOUI 6 431) (IUATA( JJ) pUJJ=1°91) 
LE (UJde¢Qe1IGO TO 445 

LE (Jd e EQ 0 4) WRITE CLOUT 2 432) (IVATA(UJJU) ep JJJR1 991) 
ah (Ud eh O24) G0 TO 445 

LE (Ud e Qe 7) WRITE (CLOUT ® 433) (IDATA (JUNI) e JJJ=1°91) 
LF (JUe=Qe7)G0 TO 445 

LF ve EQe LJ) WRITE (ILOUT #434) CLDATA (JJ) e JJUH1°91) 
ak (JdelQelLUIGU TO 445 

WKITE (CLOUT 2435) (IDATA (JU) e JJJ=1 0°91) 
rURMAT(LX*oHLUt1L Yr9iAl) 

FORMAT (LX*eonlLUm-2 Yr9uAl) 
FUORMAT(LXeotibUneS Yr 9lAl) 
FUOKMAL(LX*o0HlU-& Yr9LAl) 

FORMAT] (10H Y*9lAl) 

CONT aiNUE 

CUuNTIWwue 

WRITco (CLOUT? 437) 

FunMAl (0A°9201HIX)) 

Viki Toe (CLOUT *® 439) 
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FORMAT (7A 0 1H0 9 8X 9 QHLU e BX 2H20 6 8X9 QH3U 9 BX 9 2H4U 9 BX 9 CHSUPBXe ZHOU 
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WRITE CLOUT? 433) 

FUOKMAT(4UX*e LSHANGLE OF INCIDENCE) 

JO 491 JJ=1°13 

J=7*#JU=—0 

JJJ=7* 40d 

Wale’ 70492) (ROGL (UK) e JKK=Ue UUJ) oe Me UU 
FORMAT (2X0 7F9ebe 3Ke ler 2Xe 2) 

COnmT TIWUE 

CONT aWwUe 

AK LTcVLOUT »®250U) 

WRITE VLOUT*? 505) 

WRITE (TOUT? 510) 

UY blu JJ=199 

JIJJ=JUIU*" 10-9 

JJIIJ=IJUItY 

JIJUIJJ=JI¥* 10-10 

WL TEC SOUT® 311) JIJJJe (RTOT(LL) eLLoUJUe WU) 
COWTIWJE 

WxiTco(LOUTe S3L2)RTUT(YL) 

AKITE(,OUT» Sub) 

WRKi Te (CLOUT? 5u8) 

FORMAT (45Xe LUHSUMAARY OF) 

WkaTet,OUT*? 505) 

WRLTLEC {OUT + 507) 

FORMAT (4UX*e21HFOR ALej SUN POSITIONS) 
WRITE ( LOUT*® 430) 

LU 64+uU JJ=19°10 

YU GSU J=1°91 

JA=1Ll-UdJ 

ROvV=ERTOT (JU) *1-0c0+0'9 
LROG=1NT(L00*®ALOG1U(KOG) ) 

IF (IROGeGEeJA*100) LOATA( JU) =IBX 

IF (IRUGeLT-eJA*100) IDATA(JU)=IBL 

COWT LiWUE 

lr (UJeEUeLIWKITE (IOUT ® 431) (IVATA( UU) »p JJJ=1°9)1) 
IF (UJerQe-1)GU0 TO 645 

Ir (Udo rue 4) WRITE (LOUT 2432) (IUATA( UJI) »>JJJ=1 0°91) 
lFiddet We4) 00 TO 045 

IF (Ue Qe 7) WRITE C LOU 0433) (IUATAC QJ) e JWJJ=H1°91) 
TF (Jef Qe7)GU TO v645 

aF (Jd er Vel) WRITE CLOUT ®& 434) (iL aATA (UJI) »JJJ=1°91) 


A-L, 
UNCLASSIFIED 


UNCLASSIFIED 
cy 


Clue aF (Ider eli) G0 TO 645 

ele witiTe(lOUT® 455) (IDATA(JJJ) ep JJJ=1°91) 
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TYPICAL RESULTS 
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